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Abstract 


Computational Fluid Dynamic studies of a rigid parawing at Mach numbers from 0.8 to 
4.65 were carried out using three established inviscid, viscous and independent panel method 
codes. Pressure distributions along four chordwise sections of the wing were compared to 
experimental wind tunnel data gathered from NASA technical reports. Results show good 
prediction of the overall trends and magnitudes of the pressure distributions for the inviscid and 
viscous solvers. Pressure results for the panel method code diverge from test data at large 
angles of attack due to shock interaction phenomena. Trends in the flow behavior and their 
effect on the integrated force and moments on this type of wing are examined in detail using the 
inviscid CFD code results. 


Introduction 

The use of parawings for supersonic applications presents both appealing benefits and 
potential challenges. The particular concave shape of the wing poses distinctive aerodynamic 
behaviors, especially at supersonic conditions. In order to predict parawing aerodynamic 
performance at these conditions, it is important to develop an understanding of the ability of 
modern Computational Fluid Dynamic (CFD) tools to effectively model these behaviors. 

This study presents CFD results of a rigid parawing model at low-subsonic to high- 
supersonic Mach numbers using three established inviscid, viscous and independent panel 
method codes, and a comparison of the computational results to experimental data obtained 
from previous wind-tunnel studies [1,2,3]. The simulation codes tested consist of CBAERO, 
CART3D, and STAR-CCM+. The CBAERO software package is an engineering level tool for 
predicting the aerodynamic and aero-thermodynamic characteristics of general vehicle 
configurations. For subsonic Mach numbers, CBAERO uses a fast, unstructured, multi-pole 
panel algorithm coupled with a streamline tracing formulation to define the viscous surface 
solution. For supersonic and hypersonic Mach numbers, CBAERO uses various independent 
panel methods coupled with a streamline tracing formulation, an attachment line detection 
method, and stagnation-attachment line heating models to define the viscous aero-thermal 
environment. [4,5]. CART3D is a high-fidelity inviscid analysis package for conceptual 
preliminary aerodynamic design. The package uses an automatic generation of topologically 
unstructured, adaptively refined Cartesian meshes and an adjoint solver [5-7], Viscous results 
were also obtained using the commercial code STAR-CCM+ [8]. STAR-CCM+ solves the 
Navier-Stokes equations using a prismatic unstructured mesh near the wall surface, and an 
unstructured field mesh to discretize between the outermost layer of the prismatic mesh and the 
outer boundaries. Boundary layers are modeled using wall functions with a variety of 
conventional turbulence models. 

The geometry analyzed consists of a 61.6 degree swept parawing with a cylindrical 
leading edge and slender keel. The canopy forms half of a right circular cone and is considered 
to be rigid both in the computational and experimental studies. Pressure distributions along four 
chordwise sections of the wing were compared to experimental wind tunnel data gathered from 
NASA technical reports [1-3]. The reports provide wind tunnel data for a 14.6 inch rigid model 
of the same shape, at Mach numbers from 0.8 to 4.65 and angles of attack of 0 to 90 degrees. 

Description of Geometry 

A diagram of the model and the pressure port locations is presented in Figure 1. The 
experimental model presented in the reports consists of two 1/10 inch thick metal conical 
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canopy shapes so that the trailing edges forms half circles and the leading edges had 62.5 
degree sweep. The canopies were welded to a 14.2 inch center keel and plastic leading edges 
of 3/8 inch diameter were cast to the canopies to complete the model [1]. An auxiliary strut 
extended below the reference plane of the models as shown in the left view of Figure 1. Wires 
were attached from the auxiliary strut to the tip leading edge of the wings in order to reduce the 
deflection of the wings under aerodynamic loads (see Figure 1 b of Ref. 3). 

The computational models recreated all the features of the model, including the auxiliary 
strut, but excluding the wires, as their effect on the surface pressure in the wings was assumed 
negligible. Although previous studies have suggested that the shape of the flexible canopy can 
have a significant impact on the aerodynamics of an actual parawing at supersonic and 
hypersonic conditions [9], a rigid parawing model is appropriate for this investigation, as this 
study focused solely on the aerodynamic modeling capabilities of various CFD tools. The 
experimental setup is assumed sufficiently rigid, so any effects of aeroelastic deformations 
during the testing were ignored. The tunnel walls were not modeled in the computational study. 
Instead, a wall correction supplied by the experimental data was applied and will be discussed 
in the Methodology section of this report. 

Pressure was measured at specific locations on the surface of the wings of both 
experimental and computational models. These locations consist of four chordwise sections 
distributed uniformly along the span, containing 25 measurement locations each, distributed 
along the top and bottom of the wing as shown in Figure 1 . 

Methodology 

Tables 1, 2 and 3 present the test matrices for this study. The selection of test Mach 
number and angles of attack was driven mainly by the available experimental data. A significant 
level of experimental uncertainty exists on these numbers due to limitations in the information 
presented in the experimental reports. The Mach number for the subsonic case was calculated 
by using given approximate values of dynamic pressure and Reynolds number, and assuming a 
flow temperature of 20° C. A 0.15° angle of attack correction was applied to the subsonic 
nominal angles of attack based on a reported estimated sting deflection of around 0.1 °to 0.2°. 
The same correction was applied to the transonic cases. The experimental reports also stated 
that the Mach 1.2 data might have been affected by reflected compressions and expansions 
from the wind-tunnel walls, and the data was not corrected for these effects. Finally, for the 
supersonic cases, a flow misalignment angle was provided for each Mach number and was 
used accordingly. No bounds on measurement errors or wind tunnel condition uncertainties are 
known. 

The CART3D test matrix covers Mach numbers from 0.8 to 4.65 and angles of attack of 
0° to 90° in increments of 10°. An additional set of cases with angle of attack increment of 2.5° 
were carried out for the Mach 3.95 and 4.65 cases between 50° and 70° in order to study in 
more detail the flow behavior in this regime. These cases were solved using the CART3D 
adjoint solver with a total drag functional. Ten adaptation cycles were used, except in 80° and 
90° angle of attack cases where the number of adaptation cycles was limited to 7 due to 
convergence issues. 

The CBAERO tests focus on reproducing the available high-supersonic experimental 
data. The tangent-wedge surface pressure approximation option was used for these cases, as 
it was found to reproduce experimental results better than the Modified Newtonian model. 

The STAR-CCM+ study was focused on selected cases in the subsonic, transonic and 
supersonic regimes. A fully turbulent boundary layer was assumed to exist over the entire 
model, and the k-co SST turbulence model was used for all cases. 
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Presentation of Data 


Surface pressure data for experimental and computational results of the Parawing 
(Rogallo Wing) are presented in Figures 2 to 4. Figure 2 contains the results for the CART3D 
test cases compared to the experimental results. The center of the diagrams show top (left) and 
bottom (right) views of the wing, colored by computational pressure coefficient (Cp) magnitudes. 
Chordwise lines specify the location of the four measurement sections, and circles colored by 
experimental Cp measurements are located along these lines, corresponding to their location on 
the actual model. Plots of surface pressure coefficient versus normalized chordwise location for 
each of the four chordwise sections are located on the four corners of each diagram. The solid 
black lines represent computational results, while the red circles correspond to experimental 
measured data, where available. 

Figure 3 contains the STAR-CCM+ results compared to CART3D results and to 
experimental data. The center of the diagrams contains surface Cp magnitude contours for the 
bottom of the model with STAR-CCM+ simulations shown on the left and CART3D simulations 
on the right. The Cp plots on the corners show STAR-CCM+ results in blue, CART3D results in 
black, and experimental data in red symbols. 

Figure 4 shows all the CBAERO results compared to CART3D results and to 
experimental data. The center of the diagrams contains surface Cp contours for the bottom of 
the model with CBAERO simulations on the left and CART3D simulations on the right. The Cp 
plots on the corners show CBAERO results in blue, CART3D results in black, and experimental 
data in red symbols. 

Finally, Figures 5 and 6 show contour plots of pressure coefficient for a mid-keel section 
of the wing, and force coefficient results, respectively from CART3D simulations. 

Discussion of Results 

Comparison of CART3D Results to Experimental Results 

The CART3D results show an overall good prediction of surface pressure magnitudes 
and trends across the subsonic, transonic and supersonic regimes. The subsonic results (Mach 
0.2) show a very good comparison to experimental data, except for a marked over-prediction of 
pressure coefficients in the top surface of the wing at the angles of attack of 50° and 60°. This 
is consistent with the possibility of an inaccurate prediction of the flow separation angle due to 
inviscid assumptions. The transonic results (Mach 0.8 to 1.2) show an excellent correlation of 
the surface pressure trends, but a consistent over-prediction of surface pressures, particularly 
for the Mach 1 .2 cases. A possible explanation is the influence of reflected shocks from the 
tunnel walls, mentioned in the experimental reports. The supersonic cases (Mach 2.29 to 4.65) 
also show excellent predictions of pressure magnitudes and trends. 

Comparison of STAR-CCM+ Results to CART3D and Experimental Results 

STAR-CCM+ results proved to be very similar to CART3D results across all Mach 
number regimes. In the subsonic regime, STAR-CCM+ showed slightly better prediction of the 
flow separation angle when compared to CART3D. The transonic results show slightly different 
predictions for the pressures in the top of the wing, but also show a consistent over-prediction of 
pressures when compared to experimental results. In the supersonic regime, CART3D and 
STAR-CCM+ showed equal predictive capabilities. 
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Comparison of CBAERO Results to CART3D and Experimental Results 

The CBAERO results show a good agreement with experimental data and with CART3D 
results up to an angle of attack of around 30°. At this point, the results start diverging 
significantly as the angle of attack increases, particularly for the lower surface of the wing. At 
70° angle of attack the experimental data shows a fairly uniform pressure distribution under the 
wing, while CBAERO predicts a distribution that increases in magnitude from leading to trailing 
edge. As would be expected, these discrepancies can be attributed to the effects of shock 
interactions on the wing, which are not captured by the panel method. 

Parawing Supersonic Aerodynamic Behavior 

The CART3D solutions were used to examine the predicted flow behavior at supersonic 
conditions in more detail, and to determine how these predictions are supported by 
experimental data. As shown in Figure 5, the pressure contours on a transverse section of the 
wing at Mach 4.65 show a shock structure composed of three bow shocks, formed by the 
leading edges and the keel. As the angle of attack of the vehicle increases, the bow shocks 
grow larger and stronger and eventually merge together into a single shock covering the 
parawing. 

The effect of these shocks on the surface pressure distribution can be observed in the 
surface pressure contour plots for the Mach 4.65 case (Figures 2.97 to 2.112). At angles of 
attack of 0° to 40° the shocks are still weak and the pressure distribution is fairly uniform. 
Beyond 50° angle of attack, two contour lines of localized high pressure can be observed on 
each wing, parallel to the leading edge and the keel. As the angle of attack increases, these 
areas of high pressure move closer to the leading edge and keel. At 65° angle of attack, the 
lines of high pressure have reached the leading edge and the keel, and the pressure under the 
wing becomes relatively constant. This point corresponds to the merging of the three shock 
structures. From 65.9° to 90.0° angle of attack the pressure under the wing stays uniform while 
increasing in magnitude. The experimental data for these cases clearly shows the same 
behavior, indicating that this effect is physical. 

The lift, drag and moment coefficients computed by CART3D are shown in Figure 6. A 
change in the trend of these curves can be observed at around 65° angle of attack for the 
supersonic and hypersonic cases, particularly in the lift and moment plots, which corresponds to 
the shock structure merge phenomenon. The most notable effect is in the moment coefficient 
plot, where the curve abruptly changes from negative to positive slope. This could have 
negative implications in the stability of this wing at these angles of attack and regimes. The lift 
curve also shows a change in slope from positive to negative at this angle of attack. Comparing 
the pressure distributions of the panel method and the experimental results at 70°, it can also be 
concluded that the shock merge effect increases the pressure under the wing significantly, and 
this could produce an increased drag performance. 

For the Mach 3.95 and 4.65 cases, these conclusions can be directly supported by the 
experimental data, which shows very similar surface pressure distributions which would lead to 
similar force coefficients. Although there is no experimental data for the Mach 10 hypersonic 
case, the computational results predict that the same behavior is present, although it is less 
abrupt and happens at a higher angle of attack. 

Conclusions 

Computational Fluid Dynamic results of a rigid parawing at Mach numbers from 0.8 to 
4.65 were presented for three established inviscid, viscous and independent panel method 
codes. Results show good prediction of the overall trends and magnitudes of the pressure 
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distributions for the inviscid and viscous solvers. No particular advantages were observed for 
the viscous solver in comparison to the inviscid solver at supersonic conditions. Surface 
pressure results for the panel code diverge significantly from test data at large angles of attack 
due to the importance of shock interaction phenomena. Close study of the inviscid solver 
results showed a shock structure composed of separate bow shocks on the leading edges and 
the keel, which merge together into a single, large shock as the angle of attack increases. This 
process, which was clearly observed in both the computational and experimental data, has a 
significant effect on the surface pressure distributions. Based on the resulting force coefficients 
from the computational analysis, it is concluded that these effects can have a significant impact 
on the aerodynamic characteristics of the wing, particularly in a reduction in lift, increase in drag, 
and a decrease in static stability. Although no experimental data was available for a hypersonic 
case, a Mach 10 computational solution showed that these effects were still present, but 
occurred at a higher angle of attack. These results are a function of the concave shape of the 
wing and the supersonic flow conditions, and perhaps shed some light into the aerodynamics of 
other concave shapes in similar flows. Finally, these results show that modern CFD methods 
can be used to simulate these Parawing concepts and in turn help optimize them to achieve 
desirable aerodynamic performance. 
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Table 1. CART3D Test Matrix 
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• Cases where experimental data is available are marked in blue 
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Figure 1. Parawing geometry and pressure port locations 
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Figure 5. CART3D Mid-keel Section Contour Plots of Pressure Coefficient for 

Parawing at Mach 4.65 
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*Reference area is projected area from top view 
*Reference length is keel length 

*Moment coefficient reference point located at 60% keel length from front to back 

Figure 6. CART3D Aerodynamic Coefficient Results for Parawing at Various Mach 

Numbers 
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